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The evolution and self-reconnection of a perturbed trefoil vortex knot is simulated, then 
compared to recent experimental measurements ( Scheeler et aL|2014a ). Qualitative com¬ 
parisons using three-dimensional vorticity isosurfaces and lines, then quantitative com¬ 
parisons using the helicity. To have a single initial reconnection, as in the experiments, 
the trefoil is perturbed by 4 weak vortex rings. Initially there is a long period with de¬ 
formations similar to the experiment during which the energy, continuum helicity and 
topological self-linking number are all preserved. In the next period, once reconnection 
has clearly begun, a Reynolds number independent fraction of the initial helicity is dissi¬ 
pated in a finite time. In contrast, the experimental analysis finds that the helicity inferred 
from the trajectories of hydrogen bubbles is preserved during reconnection. Since vortices 
reconnect gradually in a classical fluid, it is suggested that the essential difference is in 
the interpretation of the reconnection timescales associated with the observed events. 
Both the time when reconnection begins, and when it ends. Supporting evidence for the 
strong numerical helicity depletion is provided by spectra, a profile and visualisations of 
the helicity that show the formation of negative helicity on the periphery of the trefoil. 
A single case with the same trajectory and circulation, but a thinner core, replicates this 
helicity depletion despite larger Sobolev norms, showing that the reconnection timescale 
is determined by the initial trajectory and circulation of the trefoil, not the initial vortic¬ 
ity. This case also shows that the very small viscosity, > 0 mathematical restrictions 
upon finite-time dissipative behavior do not apply to this range of modest viscosities. 


1. Background 

An intrinsic property of a field of twisted and linked vortices is its helicity 71, the 
volume integral of the helicity density h = u ■ u} (1.4). This integral is conserved by 


the inviscid equations and in a turbulent flow the helicity density can move through 
scales in a manner similar to the cascade of kinetic energy (1.2), the other quadratic 


inviscid invariant. However, a well-defined role for helicity in turbulent dynamics remains 


elusive (Moffatt 2014). Reasons for this include the difficulties in imposing physical space 
helical structures, in both experiments and simulations, as well as identifying analysis 
tools distinct from those used to study the energy cascade that can define what helical 
structures and dynamics are. 


Kleckner & Irvine (2013) have recently demonstrated one way to investigate helic¬ 


ity experimentally. This is to imprint helical vortex structures into a classical fluid by 
yanking SD-printed aerofoil knots, covered with hydrogen bubbles out of a water tank. 
The vortices shed by the aerofoil are then visualised by the hydrogen bubbles which are 
shed along with the vortices and align themselves into filaments along the low pressure 
vortex cores. Once the trajectories of these filaments are known, then estimates of the 
helicity can be determined from the topology of the filaments. This can be done by de¬ 
termining the linking, twist and writhe of the filaments, then multiplying each of these 
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Figure 1. Vorticity isosurface plus two closed vortex lines of the perturbed trefoil vortex at 
t = 6, not long after initialization. Its self-linking is Cs = 3, which can be split into writhe 
W = 3.15 and twist T = —0.1. 
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by their circulations to estimate the helicity (1.13), without knowing either the velocity 
and vorticity. 

The purpose of this paper is to simulate the most complicated configuration of|Kleck-| 
ner & Irvine (2013), a trefoil vortex, using the initial condition shown in figure Then 


address the following. First, make qualitative comparisons with the experimental trefoil 
visualisations in Scheeler et al. ( 2014a|b ) to demonstrate the relevance of these simula¬ 
tions to their experiment. Next, address their surprising claim that the trefoil preserves 
its helicity during reconnection. Third, discuss how the unique characteristics of the tre¬ 
foil make it an ideal tool for investigating the regularity of the Navier-Stokes equation in 
ways that other initial configurations cannot. 

Generating a trefoil vortex is challenging, both experimentally and numerically. The 
challenge is to weave, but not link, a vortex of finite diameter and fixed circulation 
into a (3,2) knot. That is, a knot with three crossings over two loops whose self-linking 
number (1.10) is Cs = W + T = 3. In a classical fluid, the resulting global helicity is 
exactly H = C^Cs, where T is the circulation of the vortex (Laing et al. 2015). The 


additional challenge for these simulations is to have only one initial reconnection, as in 
the experiments. Instead of the three simultaneous reconnections that an ideal trefoil 
generates ( Kida fc Takaoka ]|I987 ). 

The characteristics of the trefoil that make it ideal for investigating the regularity 
properties of the Navier-Stokes and Euler equations are the following: First, unlike other 
configurations such as initially anti-parallel or orthogonal vortices, it has finite energy in 
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an infinite domain. Second, it should self-reconnect due to how the two loops cross three 
times when the trefoil is projected into its direction of propagation. Third, opposing this 
tendency, because the initial helicity is nearly maximal, it can be used to investigate how 
helicity suppresses nonlinearities. Fourth, it can be simulated in a periodic box, making 
detailed Fourier analysis using Sobolev norms possible. 

The surprising result from the trefoil vortex knot experiment was that their measure 
of the helicity was preserved over the entire reconnection. Unlike their linked ring case, 
whose helicity noticeably decayed as the reconnection generated a gap between the vor¬ 
tices. The diagrams, but not necessarily the underlying bubble patterns, in |Kleckner| 
& Irvine (2013) and Scheeler et al. (2014a I suggest that there has also been complete 


reconnection to the simpler knots for the trefoil, even though its helicity has not changed. 

Do the simulations agree, or disagree, with these experimental result? The answer 
appears to be a bit of both, if one considers the possibility that the first reconnection event 
of the experimental trefoil has been misinterpreted. In the new simulations, the trefoil’s 
helicity is also preserved longer than the helicity for previous configurations using vortices 
with similar radii and circulations ( Kerr |2005a 2013a). However, once reconnection does 
start, viscous dissipation generates a finite, Reynolds number independent, change in the 
helicity. This might be the true signihcance to the trefoil’s unique dynamics. 

This paper is organised as follows. First, the equations and diagnostics used will be 
given. Then, the initialisation of the vortex illustrated in figureis described. Next, the 
similarities between the pre-reconnection evolution of this initial condition and evolution 
of the experimental trefoil are demonstrated using the isosurfaces and vortex lines in 
figures]^ {t = 24) and (t = 31). The remaining three-dimensional figures are then 
introduced, followed by the quantitative diagnostics that can tell us when reconnection 
begins (traditional vorticity norms in figure and when helicity decays (helicity and 
its partner norms in figure . These timescales underlie the comparisons that follow 
between the dissipation isosurfaces given in figure and the experimental movie that can 
be accessed through Scheeler et al. ( 2014b[ ). Spectra, a profile and helicity isosurfaces are 
then used to interpret the origins of the helicity dissipation in terms of a dual cascade, 
that is a cascade with oppositely signed helicity componets moving to opposite scales in 
both physical space and Fourier space. At the end is a discussion on how the observed 
finite-helicity depletion might be influenced by the current mathematical constraints 
upon the regularity of the Navier-Stokes equation. 


1.1. Equations and continuum diagnostics 


The governing equations for the simulations in this paper will be the incompressible 
(V • M = 0) Navier-Stokes equations in a periodic box and the numerical method will be 
a pseudo-spectral code with a very high wavenumber cut-off filter (Kerr |[20l^ ). 


du 

iH 


vAu 


+ {u ■ \7)u 


dissipation 


—Vp. 


( 1 . 1 ) 


The equations for its energy, enstrophy and helicity densities, e 
respectively, and their volume-integrated norms are: 

de 

— + {u ■ V)e = —V • (up) + vAe — iy(Vu)^ , 


= and h = u u) 

E=^ I u^dV . ( 1 . 2 ) 


e=dissipation 
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Figure 2. Isosurfaces at f = 24, shortly before reconnection begins. The vorticity isosurface 
is in blue and the helicity isosurfaces are of 0.05max(/i) in green and 0.05niin(h) in yellow 
where max(/i) = —0.22 and min(/i) = 0.085. The position of uJm = ||u;|joo, where reconnection 
is about to begin, is indicted by the red cross. 


dt 


+ (M-V)|a;p= 2a^ +uA|a;|2- , Z = J uj^dV. {1.3) 


Zp=production 


dissipation 


+ {u ■ v)/i = —uj ■ vn 


+ uA/i — 2utr(Va; • Vm^) K = j u ■ u)dV , (1.4) 


dh 
dt 

(^-transport ^-transport eh='H-dissipation 

where a; = V x m is the vorticity vector and W = p — ^ (pressure head pt = 

p+ \u^). The two quadratic, inviscid (u = 0) invariants are the energy E and helicity % 
( Moffatt 1 1969 1. Unlike the energy, helicity can be of either sign, is not Galilean invariant 
and can grow due to its viscous terms (Biferale & Kerr 19951. The enstrophy Z can 
grow due to its production term Zp. 

One can also determine spectra, the spectral transfer of energy and a variety of higher- 
order Lebesgue ||u||lp and Sobolev norms: 

II^IIlp = [y J dV\u\P^ and ||u||^, = (^J d^k\ky\u{k)\- 


1/2 


(1.5) 











Simulated Navier-Stokes trefoil reconneetion 


The vorticity diagnostics in this paper will be: 
LOm = sup |a;| = Halloo = ||u;|iL 


and Z = ||a;|j 2 = 




( 1 . 6 ) 


(1.7) 


plus the normalised enstrophy production, or velocity-derivative skewness, 

-Su = Cs^ 

where Cs is a constant relating this isotropic property to the original anisotropic hot-wire 
measurements of the velocity derivative skewness (S'„). 

The following two diagnostics, with the same dimensional scaling as the helicity, will 
be shown in figure 

La = ||M||i3 and iJi /2 = ||m||^i /2 • (1.8) 

Control of these norms has been shown to be critical to our understanding of the regu¬ 
larity of the Navier-Stokes equations ( Escauriaza et a/)]|2003 Seregin 


1.2. Vortex lines and linking numbers 

To provide comparisions to the experimental vortex lines, figures and [^include vortex 

lines Ci defined by trajectories Xj{s) that are tangent to the vorticity field uj{xj{s),t) at 
a given time. These are determined by the Matlab streamline function that solves the 
following ordinary differential equation: 

dxj{s) 


ds 




(1.9) 


over the arbitrary variable s from seed positions Xj{0). The vorticity vectors u:{xj(s)) 
are interpolated from the vorticity field on the Euleran mesh. The seeds were chosen from 
the positions around, but not necessarily at, local vorticity maxima. 

Comparisons of these vortex lines, their topological properties and associated vortic¬ 
ity isosurfaces with the experimental hydrogen bubble visualisations and analysis from 
Scheeler et al. ( 2014a|b | are the basis for the discussion in section 3.4 about reconnection 
and whether helicity is preserved or decays. 


The four topological numbers introduced by Moffatt & Ricca (1992) for generating 
helicity are: The intervortex linking numbers between distinct vortex trajectories 
Ci and Cj. The writhe Wi and twist 7) numbers of a given vortex. And their sum, the 
self-linking number: 

Cs^ = W^+%. ( 1 . 10 ) 

The Cij and Csi for closed loops will be integers (Pohl 19681, but that is not a requirement 
for either the writhe W or the twist T. 

The quantitative tool that will be used to determine the writhe, self-linking and inter¬ 
vortex linking will be a regularised Gauss linking integral. 


1 


- E 4^ 


{dxi X dxj) 


{\Xi - Xj\‘^ 


(52) 


1.5 


( 1 . 11 ) 


where <5 = 0 when f yf j and for the calculating the self-linking Cs, which can be obtained 
by choosing two parallel trajectories of the type illustrated in figureThese trajectories 
can be thought of as the edges of a vortex ribbon. The writhe Wi is defined by taking 
i = j in (1.11) (Calugareanu 1959) and calculated with S 0. The twist 7) can be 


determined from the line integral of the Frenet-Serret torsion of the vortex lines: 

/ rds, where r = —— 

ds 


% = 


B. 


( 1 . 12 ) 
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Growth for reconnection 



Figure 3. Time evolution of the scaled enstrophy 40^, normalized production Su and scal¬ 
ing maximum vorticity 0.1||aj||<3o for two cases. Q025 with v = 2.5 x 10“'^ and Q00625 with 
V = 6.25 X 10"®. 



Figure 4. Time evolution of the normalised helicity , where £ is the size of the period 

box, for 4 viscosities; v — 0.0005 to u = 0.0000625. By t ~ 72 all 4 cases have roughly the same 
decrease in helicity. Not shown is case R05 with a thinner core and u = 5 x 10“'*, whose helicity 
decay is almost identical to case Q025 with u = 2.5 x 10“'^. Inset shows normalised Ls and 
Hi /2 = ||w||jji /2 for two of the calculations. La, Jfi /2 and R are all normalised to have the units 
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N and B are the normal and bi-normal Frenet-Serret unit vectors along the trajectories. 
The focus in this paper will be on the self-linking Csi, the most robust of these numbers. 
£si will also provide a useful diagnostic during the later stages of reconnection when it 
became impossible to find trajectories that formed closed loops. 

Obtaining reliable values for the CsiW and T numbers was not straightforward using 
this numerical data due to the irregularities in the trajectories of the vortices determined 


using (1.9). Smoothing the resulting trajectories would help, but the following procedures 


were used instead. 

• For £ 5 , using two independently determined trajectories starting from near, but not 
identical, seed points was better than adding a push given by either the curvature N or 
bi-normal B vectors to the original trajectory due to noise in the variations of N and 
B along the trajectory. 

• Obtaining consistent writhe W numbers required using very small, but still finite, 5 


in ( 1 . 11 ). 


• The difficulty in calculating the twist is that three derivatives of Xj{s) are needed 
to calculate the torsion r. Reducing twist to a sum helps: % = (27r)“^E dN{s) ■ B{s). 

S 

However, values from neighbouring trajectories still varied by about 20%. Therefore, all 
the twist numbers in this paper will be the difference between the self-linking £5 and 
the writhe W. 

Once the linking, writhe and twist numbers have been determined, if the vortices have 
distinct cores, each with a well-defined circulations F,, then the helicity of the flow is 
exactly ( Laing et aL||2015 ) 

(1.13) 


H = 2^rj,£,,+r2£s, 


The factor of 2 on the Cij terms is needed because a link has two crossings in planar 
projections, instead of the single crossings one gets for a writhe or twist. 

While these trajectories are useful in identifying the initial topological changes due 
to reconnection, they are not used here for t > 31 because they very rarely formed 


closed loops. Therefore, the continuum helicity (1.4) will be the primary diagnostic for 


identifying changes in the topology during reconnection. Helicity can be applied even 
when the location and local alignment of the reconnecting vortices are difficult to identify. 
Helicity, instead of vorticity norms, will also be used to address regularity properties. 


The trefoil’s vorticity norms (1.6) are used to identify important timescales, but are not 


used for regularity questions because their growth is only modest. So modest that the 
energy dissipation e = vZ —> 0 as —)■ 0. 


2. Physical space initialisation 

The goal of the physical space initialisation is to map an analytically defined trefoil 
vortex onto an Eulerian (static) numerical mesh. The trefoil trajectories discussed in this 
paper are defined by: 


x(((>) = 

y{(t>) = 
z[4>) = 

where r{(j)) = 
and a = 


r{(j)) cos(a) 
r(0) sin(a) 

0.5cos(a) 

Tt + riacos(0) -I- asiiv{w(l) + 4 >q) 
4> + a cos{w4> + (fo)/(wrt) 


( 2 . 1 ) 


with Tt = 2, a = 0.5, w = 1.5, (fo = 0, ri = 0.25 and f) =[1 ■. Att], This weave winds itself 
about the following perturbed ring: {r{(j>) — 2 )^ -|- = 1 where r((/)) = rt+ riacos{(j>). 
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If one chooses ri = 0, one gets the traditional trefoil with a three-fold symmetry. The 
symmetric trefoil generates three simultaneous reconnection events at the three crossings 


of the trefoil loops, which is what the Gross-Pitaevski calculation in Scheeler et al. { 2014a) 


and the Navier-Stokes calculation in Kida & Takaoka (19871 find. 


The ri ^ 0 perturbation was added in the hope that with this minor adjustment, 
the trefoil would have a single initial reconnection, similar to what the evolution of the 


experimental, hydrogen-bubble trajectories show (Scheeler et al. 2014aI. This was not 


sufficient to break the three-fold symmetry, as explained below. 

Profile and direction Once the trajectory of the trefoil is established, then the 
surrounding vorticity profile and vorticity direction needs to be mapped onto the com¬ 


putational mesh. The basis for this mapping is described in Kerr (2013a). It starts by 


identifying the closest locations on each filament for every mesh point and the distances 
between these filament locations and the mesh points. These are then used to provide 
that filament’s contribution to the vorticity vector at the mesh points, a vector whose 
direction is given by the tangent to the filament at the specified location and whose mag¬ 
nitude is given by inserting the filament’s distance into that filament’s vorticity profile. 
The profile function is based upon the Rosenhead regularisation of a point vortex: 


|w|(r) = r 


(r2 + r2) 


2\2 


( 2 . 2 ) 


In addition, for every mesh point two points on the trefoil are needed. One point from 
each loop. To avoid overcounting, the space perpendicular to the central z axis is divided 
into octants. First, the octant containing the mesh point is found. Then the nearest points 
on the two loops are found from that octant, plus its two neighbouring octants. After the 
profile has been mapped onto the mesh, the field is then smoothed with a exp(— 
hyperviscous filter. 

First calculations For all of the calculations presented in this paper ri = 0.25, which 
yields a perturbed trefoil similar to the experimental initial condition. Unfortunately it 
does not generate a single initial reconnection as in the experiment and instead relaxes 
into the three-fold symmetry of an ideal trefoil before slowly diffusing into a single wavy 
vortex ring. 

Increasing the perturbation coefficient ri to values as high as 1.5 did not change this 
scenario, even though the initial perturbed structure looked nothing like the experimental 
hydrogen-bubble trefoils. 

Clearly some other type of perturbation is needed. It was pointed oui[t]that the plat¬ 
form that the trefoil model was placed upon probably generated additional independent 
vortices that would not have been identified by the hydrogen bubbles in the experiment. 
Therefore, four low intensity vortex rings propagating either in +z or —z were placed on 
the periphery of the trefoil to give it the type of external perturbation that the platform 
might have generated. 

Final initial condition The result of adding the four rings to the perturbed trefoil 
using ro = 0.25 and kf = 11.9 is given in figure]^ at t = 6, shortly after the calculation 
began. The two diagnostics are a isosurface of 0.55max(||w||oo) and two closed vortex 
filaments through the centre of the trefoil ( |1.9[ ). The five calculations discussed and 
listed in the table all use this perturbed trefoil trajectory, have a circulation of F = 0.22 
and use the same profile function, but one has a different radius. 

The circulation F, the pre-filter initial radius of rp = 0.25 and the spectral filter cutoff 
kf = 11.9 were chosen such that for cases QOS to Q00625, ojq = 1. To get a thinner core, 


t A.W. Baggaley and C.F Barenghi, private communication, 2015. 















Simulated Navier-Stokes trefoil reconneetion 


9 


Cases 

ro 

Win 

kf 

Wo 

Zo 

P 

Final Mesh 

QOS 

0.25 

1.26 

11.9 

1 

6.5 X 10"^ 

5 X lO"'^ 

512® 

Q025 

0.25 

1.26 

11.9 

1 

6.5 X 10"^ 

2.5 X 10-^ 

1024® 

Q0125 

0.25 

1.26 

11.9 

1 

6.5 X 10"^ 

1.25 X 10“'‘ 

1024® 

Q00625 

0.25 

1.26 

11.9 

1 

6.5 X 10"^ 

6.25 X 10“® 

2048® 

R05 

0.175 

2.5 

16.8 

1.9 

1.2 X 10"^ 

5 X 10"^ 

1024® 


Table 1. Parameters for the initial conditions, initial peak vorticity uo = ||aj||cx>(t = 0), initial 
enstrophy Zq, viscosity of the run and final mesh. The initial helicity for all of the calculations 
is Hit = 0) = 7.67 X 10“"^. The initialisation parameters are: The vorticity profile radius ro, the 
pre-filtering input peak vorticity tUin and the smoothing wavenumber kf in: exp(— 


case R05 used tq = 0.175 and kf = 16.8 with ojq = 1.9. The spectral filter increases the 
effective radius, defined as Ce// = (T/wo)^^^, such that for cases Q05 to Q00625 with 
ro = 0.25, Teff = 0.47 and for case R05 with ro = 0.175, re// = 0.33. 

The three lengths that can be seen in figure and are needed for comparisons with the 
experimental trefoil are: The radius of the trefoil rt = 2. The separation between the loops 
within the trefoil Sx = “2.0 = 1. And the isosurface is located at the pre-filtered profile 
radius of ro = 0.25. The highest Reynolds number simulation, with i/ = 6.25 x 10“®, is 
Rer = Y/v = 3500, 

Based upon a review of when reconnection begins in several past calculations, (Kerr 


2005b 


2013a) the relevant nonlinear timescale will not be Wq , but instead will be 


tx = 251IT = 9. 


(2.3) 


tx is the same for both initial conditions. 

For the trefoil experiment in Scheeler et al. ( 2014a|b ), the initial condition has rt « 
100mm, Sx ~ 50mm and one can take the initial profile radius to be that of the leading 
edge of the aerofoil, rg = 4mm. Much thinner than either of the numerical vortices. The 
experimental Reynolds number is Rer = 2 x 10^, an order of magnitude larger. The 
equivalent nonlinear timescale using (2.3) gives tx ~ 200ms, consistent with the times 


given for the figures for the Movie S4 in Scheeler et al. (2014b). 

Size of box. Since this initial state is generated in a periodic box, an important 
question is whether periodic boundaries might influence the trefoil’s evolution and ap¬ 
pearance. To test this possibility, case Q05 was run in a (47r)^ box. This did not change its 
appearance or any of the primary diagnostics. A related question is whether the energy 
or other norms increase as the periodic domain is made larger. If the initial condition 
were a simple vortex ring the energy would increase, but if one takes two colliding rings 
as in Lu & Doering I (20081, the energy has an upper bound as the periodic domain is 


made larger. Where does the trefoil sit? This is important for the regularity questions 
discussed in Sec. lU 

Figure 10 shows the effect of changing the size of the periodic box from (27r)^ to as 
large as (IOtt)^ upon several norms. Most importantly, the kinetic energy and Lg = ||n||L3 
do not grow significantly once the box is larger than (37r)^. This conclusion is supported 
by the energy profile in figure and the three-dimensional isosurface plot of the energy 
in figure [Tl] both of which show that the energy is almost entirely confined to the interior 
of the trefoil and does not have the type of jets above and below the knot that would be 
expected if this were a vortex ring. 
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The evolution of the trefoil reconnection will be divided into three parts. First, an 
helicity-invariant period leading up to the beginning of reconnection at t « 30. Next, a 
short period with only small changes to the helicity, but clear changes in the topology 
in the immediate vicinity of the reconnection. Then a longer, final period where these 
transformations are completed and during which there are significant and finite changes 
in the helicity. 

An overview of the changes to three-dimensional structure during these periods will 
be illustrated using the following five figures: Figure at f = 6 to illustrate the initial 
state. Figure at t = 24 to illustrate the changes just before reconnection. Figure at 
f = 31 to demonstrate how the reconnection begins. Figure at t = 36, which has the 
first clear signs of reconnection. And figure at t = 63, which shows helicity isosurfaces 
of both signs at the end of reconnection. Figure [TT] at < = 45 has been added to highlight 
how the kinetic energy is localised within the trefoil. 

To show how the trefoil core disintegrates over time and to provide a reference for 
the additional properties in these figures, all six figures include one relatively large value 
vorticity isosurface. These are supplemented with a variety of additional isosurfaces of 
vorticity, helicity, dissipation and energy, plus vortex lines. 

Figure at t = 36 has the most additional information. This is to support a proposal, 
outlined at the end of section |3.4[ that the topological changes identified as the first and 
second reconnection events in the experiments could instead be viewed as two phases of 
a single reconnection event. 

Most of the early time (< = 6,31,36) three-dimensional graphics are taken from case 
QOS with V = bx 10“"^. Vorticity isosurfaces from case Q025 with v = 2.5 x 10“^ showed 
only minimal differences in these qualitative large-scale features. Greater differences were 
evident at later times, so case Q025 is used for the t = 45 and f = 63 figures. 

After the qualitative overview using these visualisations, more quantitative diagnostics 
are discussed. 


3.1. Overall evolution 

Figure at t = 6 is used to illustrate the initial condition using a single isosurface and 
two vortex lines obtained using different seeds around the position of maximum vorticity. 
The self-linking number using (1.111 is Cs = 3. 

Figure]^ at t = 24 provides a planar view with the position of ujm = 11 I loo marked by 
a red cross along with two added low magnitude helicity isosurfaces. These isosurfaces 
are for 0.05max(/i) and 0.05min(/i), where max(h) = 0.22 and min(/i) = —0.085. Different 
perspectives show that to the right of ujm, in anticipation of the upcoming reconnection, 
the two loops are beginning to hook around one another and form a locally anti-parallel 
alignment. Also note the h < 0 yellow regions in the upper left and bottom near the 
other two loop crossings. These are forming in the directions away from red cross and 
are a result of the w-transport term in the helicity equation (1.41, with h < 0 moving 
in one direction and h > 0 moving in the opposite direction. While still weak because 
|max(h)| >> |min(/i)|, this transport of helicity along the vortex lines, h > 0 moving 
towards the dissipation sites around the reconnection and h < 0 in the opposite direction, 
is how negative helicity eventually migrates outside the original trefoil zone, as indicated 
in figure [9]at t = 63. 

Figure^ at t = 31 shows three lines in addition to its single isosurface, with their 
topological numbers given in the caption. The green trajectory originates from near, but 
not at, X, the point of uJm-, and retains the basic trefoil structure, including circumnavi¬ 
gating the central z axis twice. The point between the closest approach of its two loops is 
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Figure 5. A single vorticity isosurface plus three closed vortex lines at t — 31. The green 
line follows a remaining trefoil trajectories seeded near ujm, indicated by X. Its Cs = 4, which 
can be split into W + T = 2.85 + 1.15 = 4 The orange cross is the “reconnection point”, 
the point between the closest approach of the trefoil’s two loops and where, due to an ex¬ 
tra twist, the loops are locally anti-parallel. The Red Cs = 0 and blue Cs = 1 lines orig¬ 
inate on either side of the reconnection point and are linked, which gives a total linking of 
Lt = 2£rb + C-Sb + C-Sr = 23-1-1-0 = 3, the linking of the original trefoil. 


indicated by the orange cross, about which these two segments have adopted an anti¬ 
parallel orientation, consistent with the argument for why the total linking number might 
be preserved during anti-parallel reconnections (Laing et al. ]M^ . Now consider the line 
between these anti-parallel points and through the orange cross. The red and blue 
trajectories are seeded from points near the orange cross and in opposite directions 
perpendicular to this line. Each of these trajectories circumnavigates the central axis 
only once before returning to its origin, making them rings. Rings that stay within the 
envelope of the original trefoil and are linked through one another as they pass through 
the reconnection zone. 

On this basis, the orange cross, and not X should be viewed as the reconnection 
point. The twists of opposing polarity near X on the red and green trajectories are the 
basis for the comparisons to S4 experimental trefoil movie (Scheeler et al. |2014b[ ) in Sec. 


Figure at t = 36 uses a variety of isosurfaces to show the immediate effects of the 
reconnection. The view has been rotated by about 90° with respect to the view in figure 
[^so that the reconnection zone can be seen directly and an evolved version of the green 
trefoil from figure is included to help orient the viewer. There are two isosurfaces of 
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vorticity, the lower one at O.llujm in cyan, to show the envelope of the original trefoil, 
the second at 0.28ujm in green to emphasise where vorticity remains large. It and the 
helicity isosurface of 0.15max(/i) in blue also show a gap near [x,y] = [5.5,4.75], which 
is where the reconnection began at t = 31. The dominant large-scale vorticity structure, 
indicted by the O.llojm isosurface, is still a trefoil. 

Added to these are isosurfaces of dissipation and one of helicity advection at 0.5 max((t 6 - 
V)h) in purple. The dissipation isosurfaces are for enstrophy dissipation, 0.5max(etj) 
(|1.3[), in orange, and for helicity dissipation (1.41 of both signs, h = 0.5max(/ig) in 


red-brown and h = 0.5 min(/i) yellow. These are discussed in section 3.3 


Figure at < = 63 shows a later stage of the reconnection process. To support the 
evidence from helicity spectra and profiles for the creation of large-scale negative helicity, 
a planar view is used. The dominant vorticity structure has not been investigated in 
detail, but vortex lines suggest that what remains of the trefoil has been transformed 
into a single twisted ring. 

3.2. Quantitative diagnosties 

To complement this qualitative picture of the reconnection process, quantitative vorticity 
and helicity-related diagnostics are provided by figures [3|an d 

Figure 1^ gives three vorticity norms: ||a)||oo and Z p]. 6 [ ), which are scaled, plus the 
production skewness —Su (1.7), which is not scaled. Su is not scaled because its unsealed 


values can tell us the following: It has long been known from both experiments and 
simulations that —Su ~ 0.5 — 0.7 in fully developed turbulence. What is less appreciated 


is that for almost all initial conditions since the original pseudo-spectral paper (Orszag 


& Patterson 1972), —Su initially overshoots these values. And based upon how strong 
the overshoot is, —Su can indicate whether finite energy dissipation as the viscosity goes 
to zero is possible, or whether the dissipation goes to zero as z/ —>■ 0 . 

For the trefoil, all of these vorticity norms have significant growth over the period of 
these calculations, growth that increases as the viscosity decreases. The strong growth 
starts around t = 30, setting the timescale for when reconnection begins, and continues 
for all three norms continues until t « 50 in all cases. At that time the production 
skewness saturates at a value of —Su ~ 2.5 in the case Q00625, v = 0.0000625. Telling 
us that for t > 50, the enstrophy Z cannot grow more than exponentially fast. 

While —Su = 2.5 might seem large, compared to isotropic turbulence values of — iSu « 
0.7, the dependence of this saturation value upon the viscosity is weak. So weak that less 
than 1% of the energy is dissipated in the v = 0.0000625 calculation. That is, for th = 72 

AEu= [ dte= [ dtvZ < l%E{t = 0). 

Jo Jo 

While these vorticity norms have only modest growth, the helicity in figure [^ tells a 
different story. There are two trends. i 7 i /2 and L 3 (1.8) are nearly steady, with H 1/2 
increasing about as much as one would expect from its y/2EZ upper bound. L 3 even 
decreases slightly. Even helicity remains steady up to when reconnection begins at t « 30, 
and even longer for the higher Reynolds numbers. 

However, once reconnection has clearly formed at t « 36, the numerical helicity changes 
dramatically with H dropping by ~ 25% in all cases by th ^ 72. Unlike the conclusions 
of the experiments. 

Put another way, as the viscosity decreases in these calculations, the terms responsible 
for the viscous helicity dissipation increase. Consistent with all earlier calculations of how 
helicity changes in simulations ( Holm fc Kerr |[2007| Kimura fc Moffafr]|2014 ). 

This conclusion is supported by the manner in which an approximate non-integer self- 
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Figure 6. The perspective for f = 36 is rotated 90° clockwise from the t = 31 figure so that 
the reconnection gap to the right of the orange cross and between the green ui and helicity 
isosurfaces can be seen directly. The vorticity isosurfaces are at 0.33um cyan and 0.53cum green. 
The only closed trefoil line found is the green vortex line and is shown to help one relate this 
orientation to that at t = 31 from in figure!^ Additional isosurfaces are: blue for h = 0.15/1^; 
orange for 0.5ei.jm=0.5max(e(j) where maxle^) = 185; red-brown for 0.5ehm+=0.5max(e(,) 
and yellow for 0.5£/im-=0.5min(eji), where max(eh) = 4.9 and min(e;i) = —3.2; purple for 
0.5max((u • V)/i), helicity-transport where max((u • V)/i) = 0.06. There are black and orange 
crosses at the positions of LJrn and max(e„) respectively. 


linking number Cs decreases in time. This £5 has been calculated by applying the Gauss 
linking integral ( 1 . 11 ) to pairs of trajectories that are both seeded near the position of 
||w||oo and circumnavigate the central 2 ; axis twice, but do not close upon themselves. 
For t = 36 this gives £5 = 3.5, then 1.3 < £5 < 1.8 for 39 ^ < 48, and another jump 

down to £s < 1 for t > 66 . 

The initial structural changes in the numerical reconnection are discussed next in 
more detail, followed by possible reasons for the differences between the experimental 
and numerical conclusions in section 13.41 


3.3. Initial reconneetion period in the simulations 

To illustrate how the helicity dynamics is connected to reconnection the following ad¬ 
ditional isosurfaces have been added to figure at t = 36: A low threshold vorticity 
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isosurface, a middle-level positive helicity isosurface, three dissipation isosurfaces at half 
the maxima of the dissipation of enstrophy, maximum of the helicity dissipation, and min¬ 
imum of the helicity disspation. Plus a helicity transport isosurface at half its maximum 
and a single closed trefoil trajectory, summarised here: 

• The |a;| = O.SSwm isosurface and the trefoil trajectory show that most of the trefoil 
is still intact. 

• The helicity density h is still almost entirely positive, with the blue isosurfaces of 
h = 0.15max(/i) coincident with about half of the strong O.SSwm vorticity isosurfaces. 

• Strong enstrophy dissipation is continuing, that is reconnection is continuing within 
the orange isosurface of 0.5max(e(^) that is sandwiched between the two helicity dissi¬ 
pation isosurfaces of opposite sign just to the right of the X. 

• There is strong positive helicity dissipation in two regions, indicated by red-brown 
0.5eft+ = 0.5max(e/i) isosurfaces. 

- First, a small volume above x = 5, that is inside the gap in the cu = 0.53ujm and 
h = 0.15max(/i) isosurfaces where reconnection began at t = 31 in figure]^ 

- Second, a larger isosurface covering X. 

• There is also strong negative helicity dissipation, indicated by the yellow 0.5e/j_ = 
0.5min(e;j) isosurface of roughly the same magnitude as the positive helicity dissipation 
isosurface in red-brown. This region could, to some degree, compensate for some of 
the overall positive helicity dissipation and help explain why helicity does not change 
significantly until t > 40 in the largest Reynolds number simulations. 

• Since the strong helicity has now separated from the reconnection zone, how can 
helicity continue to decay? It can continue to decay if it is transported to where the 
dissipation is occuring, either by advection or the extra transport term along the vortex 
lines. 

- To illustrate this possibility, isosurfaces of positive helicity advection: 0.5max(M • 
V/i), are shown in purple. Negative advection regions are immediately clockwise (or to 
the left) of two y < 6 positive regions. In particular note the u • V/i surface just above 
the X that crosses between the region of strong positive helicity to its left and that of 
strong positive helicity dissipation over the X. This shows how the advection is feeding 
positive helicity into the dissipation zone. 

• This region of helicity transport and positive helicity dissipation is roughly coincident 
with the twists on two of the t = 31 trajectories in figure Twists of opposite polarity 
on the red ring, and a single twist on the t = 31 green trefoil, whose self-linking £5 = 4. 
Note that £5 = 4 is greater than the original £3 = 3 of the trefoil. Further analysis 
shows that this £s = 4 can be approximately split into the original trefoil writhe, that 
is yy « 3, and a new T « 1 twist. 

3.4. Reconnection period in the trefoil experiment 
In the discussion above it has been noted that, as in the experiments, the trefoil sim¬ 
ulations resist reconnection and, for the higher Reynolds number cases, preserve their 
helicity until t « 40, which is after reconnection has begun. However, as the reconnection 
proceeds further, the helicity in all cases drops significantly. 

It will be assumed in the next discussion that the question should not be if helicity is 
eventually lost in the experiments, but when this occurs. Is the helicity lost during the 
time period the experiment covers, just after this time period, or at much later in time? 
This is based upon realising that eventually, as the turbulence dies, energy dissipation, 
e = vZ, and the enstrophy Z, must decrease significantly. This implies that helicity must 
also decrease significantly due to the following bounds: \'H{t)\ ^ H ‘^/2 ^ \/‘^E{t)Z{t). 
Recalling that for all time E{t) ^ E(t = 0). 
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Figure 7. Helicity spectra 'Hn = 'Hfkn) and 0.25Hl^2ikn) spectra at t = 16, 30, 54 and 60 
from case Q0125. Hn < 0 values are shown in red. 0.25Hy2ikn) is included because it provides 
a guiding, but not strict, upper bound for |'H„| for all wavenumber shells. These particular times 
are shown because the Hn spectra are dominated by Hn > 0, the original sign of helicity, in 
all wavenumber shells with k < 200 (logj^Q(A:) ^ 2.3) except for these two time periods. For the 
16 ^ ^ 30 time period, the period prior to the beginning of reconnection when the trefoil is 

rearranging itself internally, Hn < 0 spans a range of intermediate wavenumbers. This could 
be due to the formation of writhe of opposite polarity to the twists, twists like those along the 
trefoil curve in Fig. The second period with Hn < 0 is represented by the spectra at t = 51 
and t = 61, near the end of the reconnection period. At this time, Hn < 0 is dominant at the 
lowest wavenumbers, that is the largest physical scales. 


There are several, overlapping reasons why this particular trefoil experiment is not 
seeing the depletion of helicity seen in the simulations, despite the similarities up to, and 
into, the beginning of reconnection. Possibilities include: 

1) The technique of identifying vortex lines from the hydrogen bubble traces is inad¬ 
equate. 

2) The principle of multiplying the lines by the square of the original circulation to 
get the helicity, which is appropriate for quantised vortices and classical vortices when 
initialised, should not be used for classical vortices with a continuous circulation once 
reconnection begins. 

3) The timescales in the experiment have been misinterpreted. 

4) The topological changes reported for the experiments have been misinterpreted. 
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possible because they have ignored how the circulation of classical vortices must reconnect 
gradually. 

Both 1) and 2) could be issues, but because their linked ring experiment gave the 
expected helicity depletion, this discussion will focus upon 3) and 4). 


Timescales. To be able to relate the reconnection timescale of the trefoil experiments 
to simulations, the trefoil simulations here and earlier work on anti-parallel and orthogo¬ 
nal reconnection, it is useful to give their reconnection times (when reconnection begins) 


in terms of their respective nonlinear timescales (2.3), and not in terms of ujq . Thus, 
anti-parallel reconnection begins at about 1.75tx ( Kerr |2013a ), orthogonal reconnection 
at 2.5tx ( Kerr |[T005b I and for the the present trefoil calculations with tx = 9, reconnec¬ 
tion begins near Z.Ztx- 


This timescale argument would suggest that the experimental and numerical trefoils 
are not in disagreement, at least in the sense of helicity conservation until after t > 


Topological changes. However, Scheeler et al. (2014a) also claim that they see helic¬ 
ity being preserved while the trefoil reconnects first into linked rings, then into a twisted 
loop. While in the simulations, helicity decay begins in all cases shortly after the trefoil 
has clearly begun to self-reconnect at t = 36. 


To understand these differences better, let us consider the time period during which 
begins in the experimental movie S4 that is accessible through IScheelei] 
\ and for which they give one frame at t = 637.7ms. This frame is from 
the period between their two proposed reconnections. The first is at t « 625ms, located 
where the twisting can be seen in the upper right in the t = 637.7ms frame and is the 
origin of separation of the trefoil into the linked rings indicated by the separate orange 
and white trajectories. The second reconnection is at t « 655ms, with a clear gap forming 
roughly where the orange and white trajectories zig-zag over one another in the middle 
of the t = 637.7ms frame and creating the final intertwined loop configuration. 


reconnection 


et al. (2014b 


However, if one compares the t = 614ms and t = 661ms frames from the movie, frames 
that represent the beginning and end of their reconnection dynamics, but do not have 
the orange trajectory, then the only significant change is the appearance of the gap from 
their second reconnection. A gap that is similar to the gap noted above for t = 36 in 
figure 

So let us consider the following alternative scenario for this sequence of events. 


• First, there is a partial reconnection at t « 625ms where some of the vortex trajec¬ 
tories are no longer following the original trefoil. Associated with this is the appearance 
of new twists along the trajectories that are not where the primary reconnection is oc¬ 
curring. Twists that are similar to those at t = 31 in figure near the X, but not the 
red X where the reconnection is actually forming. 

• Only later, and not at the same location, does the true reconnection with a gap 
become evident in the S4 movie. In S4 this starts at about t = 655ms and is similar 
to the gap that forms in the simulations by t = 36 in figure In figure the twisting 
around Wm seen in the vortex lines at t = 31, now manifests itself in the helicity and 
enstrophy dissipation terms around 

• In the simulations all of this is considered to be one reconnection. 

• This is because classical vortex reconnection is gradual. It progresses over a finite 
span of time and during this process twists on the edges of the reconnection event can 
form, twists that might look like independent reconnections. The details of this process 
will be described in another paper on anti-parallel reconnection. 
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3.5. Helieity spectra 

While the physical space analysis just presented can show the basic structure of the re¬ 
connection along with the associated transport and loss of helieity, it is limited in what 
it can tell us about how the helieity reaches the small dissipative scales. Spectra of the 
energy, the helieity and their transfer spectra can fill that gap. All of these have been plot¬ 
ted, but only the three-dimensional spectra of the helieity 77 (1.41 and of H ^/2 will 
be shown. These spectra are constructed by accumulating the Fourier transform spectra 
of 77(fc) = u{k) ■ u>{—k) and Hy^{k) = |fe||')i(A;)p into three-dimensional wavenumbers 
shells. That is for the (Stt)^ with kmin = 2/3 and kn = nkmin- 


Hn = 77(fc„) = ^ Hik) 

k„^\k\^k^+i 


and Hy 2 ikn) = ^ l^ll^(fc)P (3-1) 


The spectra will be interpreted in terms of two sets of Navier-Stokes calculations whose 
initial conditions were inspired by truncated shell, or dyadic, models (Biferale~||2003). In 


particular, how a helical spectral decomposition of the nonlinear Euler interactions (Wal- 


effe 19921 can be used to generate shell models with helicity-like invariants, including the 
popular GOY model ( Biferale &: Kerr ||1995 ). The models with the strongest interactions, 
and possible the strongest energy cascades, were based upon interactions between modes 
with oppositely signed helieity. 

Part of the objective of the Navier-Stokes calculations was to determine if these in¬ 
teractions can ever play a significant role in the dynamics. The set of calculations in 
( Biferale fc Kerr1|1995 ) used initial conditions with specified maximally helical Fourier 
modes of both signs from a single mid-range wavenumber shell. The spectra showed si¬ 
multaneous cascades that went to both large and small wavenumbers. Furthermore, when 
the higher-wavenumber cascade reached the dissipation regime and was annihilated by 
viscosity, large-scale helieity of the opposite sign was left behind. 


This effect was studied further in Holm & Kerr (2007) who looked at both helieity 


spectra and the helieity on vortices in physical space. Helieity transfer spectra showed 
spectral pulses of helieity moving to the highest wavenumbers as twists were being gen¬ 
erated on reconnecting vortex lines. Both the pulses and twists were then annihilated by 
viscosity, leaving behind helieity of the opposite sign. 

Could these observations be relevant to a trefoil? A configuration whose initial helieity 
spectrum has only one sign. These observations could be relevant if, as the helieity re¬ 
arranges itself prior to reconnection, negative Fourier modes form. These might allow 
strong oppositely-signed Fourier mode interactions to form. 

To demonstrate this possibility, figure shows three-dimensional helieity spectra Hn 
and 0.2577^^2(^n) spectra from case Q00625 for two time periods, with two times shown 

for each period. 77„ < 0 shells are in red, with the 0.2577/^2(^n) shells providing nearly 


perfect upper bounds to the absolute values |77„| of the shell helicities. A version of 
this observation that |77| ^ 0.2577^^2 has previously been noted for anti-parallel Euler 
calculations ([K err ]|2005aD . 

Figure uses spectra at t = 16 and t = 30 to represent the helieity conserving period 
of re-alignments just before reconnection, represented in physical space by figure at 
t = 24. The second period uses spectra at t = 54 and t = 60 to represent the period 
during which helieity is dissipated, represented in physical space by figure at t = 63. 
At low wavenumbers, the 0.2577/^2(^n) ^^e all the order of k~^, which is equivalent to 
energy spectra. This k~^ range is longer at the later times. 

In the pre-reconnection period, helieity is conserved, implying that if one wavenumber 
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shell develops either reduced or negative helicity, the helicity in another wavenumber shell 
must increase. Thus, at t = 16 shells 3, 4, 5 and 6 are all negative, but Q.2bHy^{kn) ~ k~^ 

only up to n = 4. At t = 30 the < 0 bands are at their Q.2bHli^{kn) upper bound, 
that is more negative, with a compensating extension of the positive k~^ to n = 6 , 
including with a slight hump. This transfer of positive helicity to high wavenumbers 
and the growth of negative helicity at intermediate wavenumbers has been confirmed by 
helicity transfer spectra. 

For the later period, at t = 54 there is a low wavenumber range where only one shell, 
n = 3, has negative helicity. Along with its neighbouring n = 2 and n = 4 shells having 
depleted positive helicity. By t = 60, there are two shells with Tin < 0. The helicity 
transfer spectra just before these times, at t = 48 and t = 54, show a transfer of Tin < 0 
from intermediate wavenumbers to the lowest wavenumbers. 

3.6. Negative helicity at large scales 

Can this cascade of negative helicity to large scales be seen in physical space? First one 
needs to identify the direction in which the negative helicity might appear, then identify 
the best viewing angles. Several profiles were taken, with the y profiles at t = 54 in figure 
[^indicating that the negative helicity region should be at y « 8.3. This is just outside 
of the envelope for 3 ^ y ^ 7.5 that contains the enstrophy. 

Now that both the helicity spectra and a profile have indicated the existence of a 
negative helicity region at large scales for t ^ 54 to 60, can a h < 0 region be found in a 
three-dimensional image? 

Figure shows that a substantial h < 0 region does exist by using the same set of 
helicity-|-vorticity isosurfaces as in figure from before reconnection. That is the same 
percentage of the maxima of w, h > 0 and h < 0, except that by t = 63 abs(min(/i)) is 
the same order of magnitude as abs(max(/i)). Thus, a significant fraction of space now 
has large h < 0. 

What does the vorticity isosurface tell us about the trefoil evolution at t = 63? By 
rotating the remaining blue vorticity isosurfaces alone, what remains is essentially a 
twisted ribbon of vorticity, consistent with the approximate self-linking number of £5 « 1 
mentioned at the end of section [3?^ Although some essense of the original trefoil does 
remain and because the helicity Ti has decayed to only 2/3rds of its original value in figure 
[^ one could argue that only one reconnection has benn completed. Since the graphics 
and Cs suggest likely that there have been two reconnections, the additional helicity is 
probably in the form of links and twists whose vorticity is below the threshold of the 
isosurface in figure 


4. Regularity 

From the point of view of the phenomenology of the physics of turbulence, the reported 
finite ATi in a finite time as v —>■ 0 trend would make complete sense. That is, the phe¬ 
nomenology that inherently assumes that there is finite energy dissipation in a finite time 
and which underlies most of the models for tnrbulent dissipation used in the geophysical 
sciences and engineering. However, from the point of view of the current mathematics of 
the Navier-Stokes equations this trend cannot continue all the way to ^' = 0. 

This is referring to results that show that for any initial condition which is regular 
under Euler evolution up to time t^, there exists a critical viscosity vq below which the 
Navier-Stokes solutions must converge to that regular Euler solution. This was originally 
formulated for periodic boxes ( Constantin||1986 1 where the definition of vq included the 
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Q025 i/=0.00025 t=54 



Figure 8. Helicity T-i, energy E and enstrophy Z profiles in t/ at t = 54, normalized by their 
maxima, for case Q025. This time is near the end of reconnection phase and has a strong low 
wavenumber negative helicity H < 0 spectral anomaly that is similar to those in hgure[^ 


size of the box £. It has since been reformulated by Masmoudi (2007) in Whole Space, 


that is for £ = oo, so that uq can be defined solely in terms of the inverse of a function 
of the time integrals of the norms (1.5) of the time-dependent Euler solution of this 
initial condition. This is the formulation that is relevant to these calculations because 
the trefoil’s volume-integrated norms are all finite as £ —t oo. 

These results should apply to these calculations so long as the Euler solutions to these 
initial conditions are regular. Such regularity under Euler is likely because by i = 24, 
just before reconnection begins and while the evolution is still inviscid, an anti-parallel 
re-alignment of the trefoil loops form to the right of the red X in figure]^ and also for the 
green trefoil trajectory in figurej^at t = 31. This alignment makes the latest evidence for 
Euler regularity from high-resolution anti-parallel calculations in |Kerr | ( |2013b[ ) relevant 
to these calculations. 

Therefore, if the viscosity were reduced far enough below the values used here, eventu¬ 
ally the trend for finite — A7£ helicity loss as the viscosity u decreases must be suppressed. 

However, this result does not apply to the range of Reynolds numbers used in this 
paper. While the relevant time integrals of the Euler Sobolev norms have not been cal¬ 
culated explicitly, case R05 with a thinner initial core and larger Sobolev norms shows 
that u > uq for the larger core cases because case R05 reproduces the convergence time 
of the other four calculations. Which would not have been possible if iz ^ i/g applied to 
any of cases Q05 to Q00625, 

For the ongoing high-resolution, higher Reynolds number, anti-parallel calculations 
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Figure 9. Isosurfaces at f = 63 as reconnection is ending. The vorticity isosurfaces are in 
blue and the helicity isosurfaces are of 0.05max(/i) in green and 0.05 mm(h) in yellow where 
max(h) — —0.41 and min(/i) = —0.34. The position where reconnection began is indicted by the 

red cross. 


that are extending those in Kerr (2013a), an explicit example will be presented for what 


happens to the i/ —>■ 0 trends when the v ^ vq threshold is crossed. And how that trend 
can be restored by increasing the size of the periodic box. 


5. Summary 

To summarise, let us return to the three periods of evolution introduced at the begin¬ 
ning. 

• First preservation. A long period during which the nearly maximal helicity im¬ 
pedes the dynamics while a locally anti-parallel alignment develops around one of the 
loop crossings, similar to the experimental observations. 

• Period at the beginning of reconnection: 

o As the reconnection begins at t = 31 in figure fragments of the original trefoil’s 
circulation are converted by reconnection into linking between new rings plus a twist 
on one of the new rings. 

o This is identified with what is called the first reconnection in the experimental 
movie, but for the simulations is interpreted as a region of twisting just outside the 
zone where the actual reconnection is beginning to form. 
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Boxsiz= (. Squares for thin trefoil R05 

Figure 10. Helicity H, energy E, L 3 and ||t<j||oo as functions of the initial box size £, 

normalised using V = £^. Plus squares for the t = 0 values for case R05 with a thinner core. 
This figure shows that for the trefoil, all of these norms have finite values in Whole Space. 


o In the simulations, comparisons between the frames at t = 31 and t = 36 suggest 
that there is period of re-alignment that could be associated with the beginning of 
reconnection from the original trefoil into a new topology dominated by linked rings, 
o However, signihcant changes in the global helicity do not appear until after a true 
gap between reconnected vorticity isosurfaces starts to appear at t = 36. 
o Which suggests that what is called the second reconnection in the experiment, 
actually represents the beginning of the first true reconnection, and only after that 
does the helicity start to decay. Which would have been after the last times shown 
in the experimental analysis. 

Later period 

o There is finite helicity dissipation in a finite time th « 72 as the viscosity goes to 
zero, th appears to depend only upon the properties of the initial trajectory of the 
trefoil, the separation between of the two loops, and the trefoil’s circulation T. 
o Along with helicity dissipation, creation of negative helicity at the largest scales 
is shown in figures |7|8|9| using spectra, profiles and visualisation, 
o One way to look at this is that a closed vortex loop can generate large-scale 
writhe and small-scale twist, whose self-linking sum is conserved (1.10). This could 
provide a mechanism for transferring negative writhe-helicity to large scales and 
positive twist-helicity to the small scales. After which, the twist can break off and 
be dissipated, leaving behind its large-scale negative writhe-helicity. 













22 


R. M. Kerr 



Figure 11. Vorticity, energy and Felicity isosurfaces of at t = 45. Vorticity is indicated by 
the blue, energy by magenta, large positive helicity by green and large-magnitude negative 
helicity by orange. The main point is that large values of the energy are entirely confined within 
the trefoil, without any signs of the type of jets above and below the knot that would be found 
for a vortex ring. 


o Over this period the energy dissipation goes to zero as viscosity goes to zero. This 
is associated with only modest growth in the traditional vorticity norms, 
o It is pointed out that if only one configuration is considered, then eventually this 
u —0 trend must be suppressed for any v < vq, where vq is a function of the Euler 
norms. 

o Despite this observation, this small viscosity can be decreased by using thinner 
initial vortices with larger Euler norms. Physically, if a given v is too small for viscous 
reconnection to relieve the small-scale twisting, then by using a thinner vortex, the 
twisting could continue. This will the subject of later work. 


6. Conclusion 

The question of a dynamical role for helicity in vortex dynamics has been considered in 
light of recent experimental claims that helicity is preserved during the self-reconnection 
of trefoil vortex. Specifically, the claim is that when lines of hydrogen bubbles marking 
the experimental trefoil vortex knot self-reconnect that the helicity determined from the 
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topological numbers (1.131 is preserved. And from this they concluded that the helicity 
is also preserved. 

Once a perturbed initial condition for the numerical trefoil vortex knot was generated, 
one that was not subject to internal instabilities and did not have symmetric recon¬ 
nections, both qualitative and quantitative comparisons with the experiments became 
possible. The first result was that as in the experiments, the helicity was preserved for 
a surprisingly long time. Long meaning up to and into the beginning of the first re¬ 
connection and up to 50% longer than the time for the reconnection of anti-parallel or 
orthogonal vortices to begin. 

However, once reconnection in the simulations begins in earnest, neither the helicity 
nor the self-linking topological number are preserved. It is argued that the experimental 


conclusion that the trefoil’s helicity is preserved through reconnection (Scheeler et al. 


2014a I is based upon how they identify the observed topological changes as coming from 


two reconnections. Based upon comparisons with these simulations, these two events are 
re-interpreted as different phases during the beginning of the first, single reconnection 
event, and before the helicity begins to decay. This suggests that the true significance of 
the trefoil experiment, if run longer, would not only be in how it suppresses reconnection 
and the depletion of helicity, but in how it should also see finite helicity dissipation by 
a fixed rescaled time as the Reynolds number increases. Even when there is almost no 
energy decay. 
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